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ABSTRACT 

Radiation pressure from the absorption and scattering of starlight by dust grains may be an important feed- 
back mechanism in regulating star-forming galaxies. We compile data from the literature on star clusters, 
star-forming subregions, normal star-forming galaxies, and starbursts to assess the importance of radiation 
pressure on dust as a feedback mechanism, by comparing the luminosity and flux of these systems to their dust 
Eddington limit. This exercise motivates a novel interpretation of the Schmidt Law, the Lir-L' co correlation, 
and the Lm-L' HCN correlation. In particular, the linear Lir-L' hcn correlation is a natural prediction of radiation 
pressure regulated star formation. Overall, we find that the Eddington limit sets a hard upper bound to the lumi- 
nosity of any star-forming region. Importantly, however, many normal star-forming galaxies have luminosities 
significantly below the Eddington limit. We explore several explanations for this discrepancy, especially the 
role of "intermittency" in normal spirals — the tendency for only a small number of subregions within a galaxy 
to be actively forming stars at any moment because of the time-dependence of the feedback process and the 
luminosity evolution of the stellar population. If radiation pressure regulates star formation in dense gas, then 
the gas depletion timescale is 6 Myr, in good agreement with observations of the densest starbursts. Finally, we 
highlight the importance of observational uncertainties — namely, the dust-to-gas ratio and the CO-to-H2 and 
HCN-to-H2 conversion factors — that must be understood before a definitive assessment of radiation pressure 
as a feedback mechanism in star-forming galaxies. 

Subject headings: galaxies: general, evolution, ISM, stellar content, starburst — stars: formation 



1. INTRODUCTION 

Understanding global star formation is crucial in under- 
standing galaxy evolution and the assembly of the z = stel- 
lar population over cosmic time. Observations indicate that 
only a few percent of the available gas reservoir in galaxies is 
converted into stars per local free-fall time (Kennicutt 1998; 
Krumholz & Tan 2007). In addition, models of the interstellar 
medium (ISM) suggest that energy and momentum injected 
by massive stars could act as a feedback loop by driving su- 
personic turbulence, which would cause most of the gas to be 
insufficiently dense to collapse, rendering star formation inef- 
ficient (Krumholz & McKee 2005). However, the interaction 
between star formation and the ISM is not well understood, 
and a mechanism for the regulation of star formation across 
the large dynamic range of star-forming environments has not 
yet been conclusively identified. Proposed mechanisms in- 
clude supernova explosions, expanding HII regions, stellar 
winds, cosmic rays, magnetic fields, and radiation pressure 
on dust (McKee & Ostriker 1977; Matzner 2002; Cunning- 
ham 2008; Chevalier & Fransson 1984; Socrates et al. 2008; 
Kim 2003; Scoville et al. 2001; Scoville 2003; Thomp- 
son, Quataert, & Murray 2005, hereafter TQM; Krumholz & 
Matzner 2009, hereafter KM09; Murray, Quataert, & Thomp- 
son 2010, hereafter MQT; Draine 2010; Hopkins et al. 2010). 4 

In the case of radiation pressure on dust, UV and optical 
radiation from OB stars is absorbed and scattered by dust 
grains and subsequently re-radiated as IR radiation. The dust 
grains are coupled to the gas of the ISM through collisions 
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4 ISM turbulence driven by non-stellar processes, such as disk instabilities, 
has also been proposed (Sellwood & Balbus 1999; Wada et al. 2002; Piontek 
& Ostriker 2004, 2007). 



and magnetic fields, so radiation pressure on the dust ex- 
erts a force on the gas as well (O'dell et al. 1967; Ferrara 
1993; Laor & Draine 1993; Murray, Quataert, & Thomp- 
son 2005). On galaxy scales, TQM showed that radiation 
pressure could constitute the majority of the vertical pressure 
support in dense starburst galaxies such as ultra-luminous in- 
frared galaxies (ULIRGs). Likewise, models of giant molec- 
ular cloud (GMC) disruption predict that radiation pressure is 
the dominant feedback mechanism regulating star formation 
in the birth of massive star clusters (MQT; KM09). In this pic- 
ture, gas in a marginally-stable galactic disk collapses to form 
a GMC and a central compact star cluster. When the stellar 
mass and luminosity of the cluster exceed the Eddington limit 
for dust, the overlying gas reservoir is expelled. Thus, the fi- 
nal stellar mass in individual star clusters is regulated by the 
dust Eddington limit. The centers of ULIRGs and GMC cores 
are optically thick to both UV and the re-radiated IR photons, 
which make them ideal candidates for radiation pressure sup- 
port since essentially all of the momentum from the starlight 
is efficiently transferred to the gas. Recent observations in- 
dicate that the most luminous GMCs in the Milky Way are 
disrupted by radiation pressure (Murray 2010). 

In this paper, we critically assess the theory of radiation 
pressure regulated star formation by comparing the picture 
developed by TQM, MQT, and KM09 with the available ob- 
servations of star-forming galaxies ranging from dense indi- 
vidual star clusters and GMCs, to normal spiral galaxies and 
starbursts. In $2] we describe the current model of radiation 
pressure feedback. We emphasize the deviations from the 
simplest version of the dust Eddington limit in assessing radi- 
ation pressure regulated feedback that arise from ambiguities 
in the value of the flux-mean dust opacity, and the tendency 
for low-density galaxies to have highly intermittent knots or 
hotspots of star formation across their disks. In $3] we com- 
pare data from the literature to models of radiative feedback. 
In |2] we discuss our conclusions, the major observational 



2 



and theoretical uncertainties in our analysis, and the implica- 
tions of our results. 

2. THEORETICAL ELEMENTS 

The statement that radiation pressure may be an important 
feedback mechanism in galaxies is equivalent to the statement 
that galaxies as a whole or the star-forming subregions within 
them approach or exceed the Eddington limit for dust, 

^Edd = , (1) 

where Fem is the Eddington flux, E is the surface density of 
the dominant component of gravitational potential in the star- 
forming region, and «f is the flux-mean opacity. The overall 
picture is that star-forming regions meet the Eddington limit 
and self-regulate in analogy with an individual massive star 
(TQM; see also Scoville et al. 2001; Scoville 2003; MQT; 
KM09). We would thus naively expect to test the theory of 
radiation pressure regulated star formation by taking the ratio 
of the observed flux (Fobs) to fkdd- However, a direct compar- 
ison between the simple theoretical expectation 

F oh jF Edd -> 1 (2) 

and the observations is complicated by both theoretical and 
observational uncertainties. For example, although gas is ex- 
pected to be the dominant mass component in and around 
massive star clusters in formation, it is unclear how best to 
estimate E in equation (Q~|i for unresolved galaxies or unre- 
solved star-forming subregions. Below, we consider both CO 
and HCN emission (see Figures[T]&|2]i, but the conversion 
from the luminosity in either of these molecular gas tracers to 
gas mass where the stars are forming, is highly uncertain. An- 
other uncertainty is the coupling of the radiation field to the 
gas, which is complicated due to both the non-gray nature of 
the dust opacity and the dumpiness of the gas on all scales 
(see jj2.1| below). Finally, there is an additional complication 
not readily apparent from the time-independent statement of 
equation the star formation rate (SFR) across the face of 
a large spiral galaxy is highly intermittent so that only a small 
number of subregions are bright at any time. As discussed by 
MQT and in detail below ( §12. 21 >. this intermittency can cause 
normal star-forming galaxies to be appear significantly sub- 
Eddington (f bs/^Edd C 1) when only their average properties 
are considered, but much closer to Eddington when a model 
is used to take this effect into account. 

2.1. The Radiation Pressure Force 

The coupling between radiation and gas in star-forming en- 
vironments is complex primarily because the flux-mean opac- 
ity Kp in equation ([TJ has a full range of more than 3 dex, 
depending on whether the SED of the system considered is 
dominated by UV or FIR light. However, there are two dis- 
tinct regimes: optically thick to UV but thin to the re-radiated 
FIR and optically thick to FIR. We call these the "single- 
scattering" and "optically thick" limits, respectively. 

2. 1 . 1 . Single-scattering Limit 

Regions in the single-scattering limit are optically thick to 
the UV but optically thin to the FIR (tfir ~ kfir E g /2). This 
limit applies over a wide range in surface density: 

E g < 5000 Mqpc" 2 ^ 1 /^ 150, (3) 
where kfir = ^2/dg, 150 is the Rosseland-mean dust opacity 
with K2 = k/(2 cm 2 g -1 ) (see $2.1.21 and /d g , 150 = /d g x 150 is 



the dust-to-gas ratio. In the single-scattering limit, UV pho- 
tons are absorbed once and then re-radiated as FIR photons, 
which free-stream out of the medium. 5 Since the column- 
averaged flux-mean optical depth in this limit is always equal 
to unity, the flux-mean opacity for the single-scattering limit 
is ~2/E g . The Eddington flux is then 
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10 M pc" 2 

where / gas = E g /E tot is the gas fraction and E tot = E g + 0.1E* 
(Wong & Blitz 2002). The wide range of column densities 
over which this limit is applicable implies that the average 
medium of most star-forming galaxies, some starbursts, and 
the GMCs that constitute them is single-scattering. 

2.1.2. Optically Thick Limit 

Dense starbursts and GMCs can reach the high gas sur- 
face densities E g > 5000 M Q pc~ 2 k^ 1 150 necessary to be- 
come optically thick to FIR photons (tfir > 1). In this case 
^rad ~ tfirF/c, and kfir depends on temperature (Bell & Lin 
1994; Semenov et al. 2003). The functional form of k fir nat- 
urally leads to two regimes: "warm" (T < 200 K) and "hot" 
(200 K < T < T sub , where r sub ~ 1500 K is the dust sublima- 
tion temperature). For typical numbers, the central tempera- 
ture of a massive, compact star cluster is 

k fir E F K F i R Mg M*W 
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eff 
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(5) 

where T e ff is the effective temperature, Kio = k/(10 cm 2 g -1 ), 
^3000 = 3000 ergs s" 1 g" 1 is the light-to-mass ratio of a zero 
age main sequence stellar population, M g .6 =M g /(10 6 Mq), 
andM*, 5 =M,/(lO 5 M ). 

Warm Starbursts — For T < 200 K, the Rosseland mean 
opacity increases as Kpm(T) ps k T 2 , where 
10~ 4 cm 2 g" 1 K" 2 /d g , i5o- In this case, 
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lO^Lokpc" 2 /^ 1 / 2 ^^. (6) 



Remarkably, the flux necessary to support the medium is in- 
dependent of E (TQM). 

Hot Starbursts — Intense, compact starbursts may have cen- 
tral temperatures greater than 200 K. The corresponding 
opacity is roughly constant with temperature: kfir^) ~ 5- 
10 cm 2 g" 1 /d g , 150 for temperatures 200 K < T < r suD . For 
typical numbers, 

E„ 



c thick 
Edd ' 



10° L kpc" 



r l/2 r i 
J gas J dg 



ISO- 



(7) 



,1O 6 M pc" 2 

The high surface densities necessary to enter this regime may 
only be attained in the pc-scale star formation thought to at- 
tend the fueling of bright AGN (Sirko & Goodman 2003; 
TQM; Levin 2007). 

5 Galaxies with surface densities less than ~5 Mq per 2 will be optically 
thin with respect to dust. Below this limit, the ionization of neutral hydro- 
gen will becomes the dominant source of opacity. The large cross-section 
(<thi ~ 6.3 X 10~ 18 cm 2 per H atom) implies an incredibly small surface den- 
sity (Eg > 10~ 3 Mq pc -2 ) is required for the medium to be optically thick to 
ionizing photons. These ionizing photons transfer momentum directly to the 
gas on the same order as the momentum transfer due to the single-scattering 
limit for dust. Thus, we encompass this limit and the single-scattering limit 
for dust under the same heading. 
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2.2. GMC Evolution & Intermittency 

In order to gauge the importance of GMC evolution and in- 
termittency, we adopt the simple picture presented by MQT 
that marginally stable (Q ps 1) disks fragment into sub-units 
on the gas disk scale height (h) to form GMCs. An individual 
star cluster is born, reaches the critical Eddington luminos- 
ity threshold, and then expels the overlying gas. Importantly, 
the timescale for collapse and expansion of the GMC is the 
disk dynamical timescale, f,jy n , which can be much longer than 
the characteristic timescale for the stellar population to de- 
crease in total luminosity, the main-sequence lifetime of mas- 
sive stars, fMS ~4x 10 6 yr. In this picture, a low-density star- 
forming galaxy with radius r should have ^(r/h) 2 sub-units, 
but only a small fraction £ (the "intermittency factor") should 
be bright at any one time. If each subregion reaches the Ed- 
dington luminosity for a time t M $ and is then dark, and then if 
a large number of subregions are averaged, one expects 



^ _ Mm iobs ?MS 

Mot ^Edd 2fdyn + fMS 



(8) 



where 



3n(2h) 
32GE to t 



1/2 
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/iioo = /i/(100pc), M on is the number of sub-units that are 
"on," and Mot is the total number of sub-units. For ex- 
ample, the normal star-forming galaxy M51 has a observed 
bolometric luminosity (L 0DS = 0.2L£dd) that is a factor of ^4 
larger than its intermittency-corrected Eddington luminosity 
C^Edd = 0-05LEdd)- Although the approximation that the stellar 
population is bright for a time Ims an d then dark is crude, the 
parameter £ gives us a way to judge the importance of inter- 
mittency in normal star-forming galaxies. 

Note that for higher densities, fdyn decreases and £ — > 1 at a 
critical surface density 



^ci-it ' 



3tt(2/!) 



32G(0.5f MS ) 2 



3x 10 3 M Q pc~ 2 h 



101). 
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which corresponds with a critical midplane pressure f cr i t ~ 
5 x 10~ 8 ergs cm" 3 /i 2 00 (see eq. ITSb . For E to t > E cr i t , massive 
stars live longer than the time required to disrupt the parent 
sub-unit (fMS > fdyn). MQT argue that in this regime the mas- 
sive stars continue to drive turbulence in the gas and maintain 
hydrostatic equilibrium in a statistical sense until fj^s, when 
the process then repeats until gas exhaustion. 

Several additional elements of GMC evolution are impor- 
tant in judging whether or not the star formation of galaxies 
is regulated by radiation pressure on dust. First, the GMCs 
collapse from regions of size h 2 , with total mass E g 7r/z 2 , and 
to a size Rqmc = h/<j>. This implies that the surface density of 
individual GMCs is Egmc ~ 2 S g , where <f> can be ~2-5 in 
the Galaxy (MQT). For example, taking cgmc = M^/Mgmc 
(oc Eg in the single-scattering limit), \P3000, and assuming 
KpiR ~ k T 2 (appropriate for T < 200 K), one finds that the 
required gas surface density for a GMC to be optically thick 

to the FIR is E^ 1 - 7000e G ^c,o.i M o P c ~ 2 - This GMC S as 
surface density would correspond to an average gas surface 
density for the galaxy that is (j) 2 times smaller (E g ~ 350- 
1800M Q pc~ 2 ). Thus, the medium surrounding the central 



star cluster may be optically thick to the FIR even if the av- 
erage gas surface density of the galaxy is less than the naive 
estimate given in 32.1.21 

At surface densities in excess of tfir = 1 for the GMCs, the 
models of MQT rely on the fact that the medium is in fact 
optically-thick to FIR radiation. This is in sharp contrast to 
the work of KM09, where they argue that the effective opti- 
cal depth is always ~1 because instabilities allow the radia- 
tion to leak out of otherwise optically-thick media. MQT ar- 
gue that the effective optical depth must be much larger than 
unity in the GMCs of dense starbursts and ULIRGs for ra- 
diation pressure to be effective as a feedback mechanism. In 
addition, they show that the effective momentum coupling be- 
tween the radiation and the gas can exceed the naive estimate 
based on a disk-averaged tfir by a factor of a few in systems 
as dense as the putative GMCs in Arp 220 because of the time- 
dependence of the GMC disruption process (see 34.31 ). 

3. RESULTS 

We compile data of super star clusters, normal star-forming 
galaxies, local starburst galaxies, ULIRGs, sub-millimeter 
galaxies (SMGs), hyper luminous infrared galaxies, and cir- 
cumnuclear starbursts to assess feedback from radiation pres- 
sure. Below, we test the hypothesis that radiation pressure 
is dynamically important in galaxies and star-forming subre- 
gions by comparing data to the Eddington limit (^2]) on a va- 
riety of physical scales ranging from globally-averaged prop- 
erties of galaxies to individual star-forming subregions within 
galaxies. 

3.1. IR Luminosity vs. Molecular Line Luminosity 

We show the total IR luminosity L\r as a function of molec- 
ular line luminosity 6 L' co (Figure [TJ and L' HCN (Figure |2]( for 
our sample of star-forming galaxies. 7 Lir is known to trace the 
total light from massive stars (e.g., Kennicutt 1998), 8 whereas 
L' CQ and L' HCN provide a measure of the total gas mass and 
dense gas mass, respectively. Under the assumption that the 
total gravitational potential is dominated by the gas on the 
physical scales where the stars are forming, the Eddington lu- 
minosity is related to L' co by 

4ttGc , 

^Edd = XcO L co , (11) 

K 

where Xqo is the L' co -Xo-M^ conversion factor and k is 
the appropriate flux-mean or Rosseland-mean opacity (either 

6 We used the J = 1-0 line unless only higher order lines were available. 

7 Aravena et al. 2008; Becklin et al. 1980; Beelen et al. 2006; Benford 
et al. 1999; Capak et al. 2008; Carilli et al. 2005; Casoli et al. 1989; Chap- 
man et al. 2005; Chung et al. 2009; Combes et al. 2010; Coppin et al. 2009; 
Coppin et al. 2010; Daddi et al. 2007; Daddi et al. 2009; Daddi et al. 2010a; 
Downes & Solomon 1998; Gao et al. 2007; Gao & Solomon 1999; Gao & 
Solomon 2004a,b; Genzel et al. 2003; Gracia-Carpio et al. 2008; Greve et 
al. 2005; Greve et al. 2006; Isaak et al. 2004; Kim & Sanders 1998; Knud- 
sen et al. 2007; Mauersberger et al. 1996; Mirabel et al. 1990; Momjian et 
al. 2007; Murphy et al. 2001; Neri et al. 2003; Riechers et al. 2006; Riech- 
ers et al. 2007; Riechers et al. 2008; Sajina et al. 2008; Sakamoto et al. 2008; 
Sanders et al. 1991; Schinnerer et al. 2006; Schinnerer et al. 2007; Schinnerer 
et al. 2008; Smith & Harvey 1996; Solomon et al. 1997; Solomon & Vanden 
Bout 2005; Walter et al. 2003; Walter et al. 2009; WeiB et al. 2001; Yan et 
al. 2010; Young & Scoville 1982; Yun et al. 2001. 

8 In reality, for normal galaxies a fraction of the UV and optical light es- 
capes before being reprocessed by dust, and a fraction of the IR is diffuse and 
likely not associated with star formation (e.g., Kennicutt et al. 2010, Calzetti 
et al. 2010). The UV and IR luminosities are roughly equal at a bolometric 
luminosity of Lbol ~ 10 10 Lq, but the UV luminosity is an order of magnitude 
larger than the IR luminosity at L bol ~ 1O 8 ' 5 L (Martin et al. 2005). Thus, 
we expect that galaxies with Lir < 10 10 Lq to move closer to the Eddington 
limit in 
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FIG. 1 . — IR luminosity as a function of CO line luminosity. The different symbols correspond to different rotational transitions of CO: solid circles (/ = 1-0), 
crosses (/ =2-1), squares (J = 3-2), and triangles (J = 4—3). The data a re the same in both panels. The lines in the left panel correspond to the single-scattering 
Eddington limit (solid line; assuming h = lOOpc and r = lOkpc; 32.1.U and the single-scattering Eddington limit accounting for intermittency (dot-dashed line; 
eq. [8). The lines in the right panel show the optically thick Eddington limit for our preferred value of the Rosseland-mean opacity (shaded region; kfir = 5- 
10 cm 2 g _1 /dg,i5o) and for an enhanced dust-to-gas ratio (dashed line; kfir = 30 cm 2 g _1 /dg,5o)- Note that no galaxies are significantly super- or sub-Eddington. 
We emphasize that it is not possible to determine which limit is applicable without knowing a surface density, so dense star-forming regions can be optically 
thick at L' c0 < 10' K km s pc 2 and approach the optically thick Eddington limit (see j|2.U . The single-scattering Eddington limit was calculated by adopting 
the Galactic CO-to-H 2 conversion factor = 4.4 M Q (K km s~ 

r 



pc ) . The optically thick Eddington limit was calculated by adopting the ULIRG CO-to-H2 



= O.8M0(Kkms-' pc 2 '"" 1 



conversion factor X^q 1rg 

single-scattering or optically thick; see jj2.1.Uj2~1.2l i. Al- 
though counterintuitive, the single-scattering Eddington lumi- 
nosity lies below the optically thick Eddington limit because 
the dust opacity is column-averaged and highly non-grey (see 
CO- We adopt = 4.4 M (K km s" 1 pc 2 )" 1 for normal 
galaxies (including a correction factor of 1 .36 to account for 
He; Strong & Mattox 1996; Dame et al. 2001), and X^ ma = 
0.8 M Q (K km s" 1 pc 2 )" 1 for galaxies with L m > 1O U L 
(as appropriate for starbursts and ULIRGs; e.g., Downes & 
Solomon 1998). Similarly, if the majority of the IR luminosity 
comes from regions where the dense molecular gas dominates 
the potential, then Lead is related to L' HCN by 



AttGc , 

^Edd - ^HCN ^HCN ) 

kfir 



(12) 



where we explicitly write n = kfir because the critical den- 
sity for HCN emitting gas is large enough that these re- 
gions should always be optically-thick ( {32.1.21 i. 9 In eq. Q~2] 
we take an L^^-to-M^™ 56 conversion factor of Xhcn = 
3 M Q (K km s _1 pc 2 ) -1 , but we caution that X^cn is uncertain 
to a factor of ~3 (Gao & Solomon 2004a,b; see j4~5l 

In Figure Q] the lines indicate the Eddington limit for var- 
ious limiting cases. The lines in the left panel show the 
single-scattering Eddington limit (solid line) and the single- 
scattering Eddington limit accounting for intermittency (dot- 
dashed line) assuming h = 100 pc and r = lOkpc (eq.fTTIi. The 
shaded region in the right panel is the optically thick Edding- 
ton limit (eq. QT|with k fir = 5-10 cm g _1 ) and the dashed 

9 Using p cr itHCN ~ 10~ 19 g cm" 3 , tfir ~ KFiRPcrit,HCN# is larger than 
unity for scales R > pc and kfir ^ few cm 2 g~' . 



line shows the optically thick Eddington limit for an enhanced 
opacity (k fir = 30 cm 2 g _I / dg)50 ; where / dgi50 = / dg x 50) due 
to an assumed higher dust-to-gas ratio in dense star-forming 
environments. We plot the single-scattering (left panel) and 
optically thick (right panel) Eddington limits separately for 
clarity, but the data are the same in both panels. The different 
symbols indicate various rotational transitions of CO: solid 
circles {J = 1-0), crosses (J = 2-1), squares {J = 3-2), and 
triangles (J = 4-3). 

Note that no galaxies exceed the optically thick Edding- 
ton limit and most galaxies are neither significantly super- or 
sub-Eddington with respect to the single-scattering Eddington 
limit. We caution that the applicable Eddington limit for any 
individual galaxy cannot be determined in this plot due to the 
lack of surface density measurements, which dictate the op- 
tical depth to the FIR and the relevant Eddington limit. For 
example, high surface density star-forming regions can be op- 
tically thick at L' co < 10 9 K km s" 1 pc 2 and lie to the left of the 
single-scattering Eddington limit (solid line in left panel) but 
below the optically thick Eddington limit (shaded region in 
right panel). For the single-scattering limit, our assumption of 
r = 10 kpc is accurate to a factor of ^5 for most of the sample, 
but the single-scattering opacity scales as r" 2 , so it is only ac- 
curate to a factor of ~25. Some compact starbursts have radii 
much smaller than our assumed radius, so they are optically 
thick, even at low L' co and this explains why they exceed the 
single-scattering limit in the left panel but are below the opti- 
cally thick limit in the the right panel. In addition, note that 
the optically-thick Eddington limit is a hard upper bound to 
a galaxy's IR luminosity, which suggests that radiation pres- 
sure feedback may set the maximum SFR of a galaxy. In the 
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FIG . 2 . — IR luminosity as a function of HCN line luminosity. The different 
symbols correspond to different rotational transitions of HCN: solid circles 
(/ = 1-0) and crosses (J = 2-1). We show the optically thick Eddington 
limit for our preferred value of the Rosseland-mean opacity (shaded region; 
"FIR = 5-10 cm 2 g -1 /dg,i5o)- The dashed line shows the effect of a factor 
of 3 increase in the dust-to-gas ratio for the optically thick Eddington limit 
(«TIR = 30 cm 2 g /d g .5o)- Note that all galaxies are within ~1 dex of the 
optically thick Eddington limit, which suggests that radiation pressure may 
regulate star formation. The optically thick Eddington limit was calculated by 
adopting the HCN-to-H2 conversion factor Xhcn = 3 M0(Kkm s _1 pc 2 ) -1 . 

Eddington-limited model, the scatter in the L iR-L' cn relation 
may be due to variations in h, Xqo (see jj4.51 >. the dust-to-gas 
ratio/metallicity (see j|4.4t . the effective radii, and the depth 
of the stellar potential. 10 

The intermittency of star formation will likely affect the 
Eddington limit for CO-emitting gas (dot-dashed line in left 
panel). L' co traces the total molecular gas reservoir includ- 
ing the molecular gas that is not actively participating in 
star formation, such as GMC envelopes and diffuse inter- 
cloud gas. The gas mass relevant for the Eddington limit 
may be overestimated for galaxies in the single-scattering 
limit. To account for this, we multiply the Eddington lu- 
minosity by the intermittency factor for CO-emitting gas 
£ - 0.06 for the Milky Way value of X co , h = 100 pc, r = 
10 kpc, and L' co = 10 9 K km s" 1 pc 2 (see eq. 1 and j2~2l . 
The intermittency factor approaches unity when L' co ~2x 
10 11 Kkm s" 1 pc 2 /zioo^io/^ccT- Thus, compact star-forming 
regions, such as the nuclear starbursts of ULIRGs, have £ ~ 1 
at low L' co due to their very small radii (e.g., Downes & 
Solomon 1998). 

Figure[2]shows the Lir-L hcn relation for our sample of star- 
forming galaxies. The shaded region represents the optically- 
thick Eddington limit (eq. [12] with kfir = 5-10 cm 2 g" 1 ) and 
the dashed line shows the optically thick Eddington limit for 
an enhanced opacity (kfir = 30 cm 2 g" 1 /d g ,5o; where /dg.50 = 
/dg x 50), which may result from a higher dust-to-gas ratio 

10 Note that previous work by Krumholz & Thompson (2007) and 
Narayanan et al. (2008) explains the slopes of the Lir-L^q and Lir-L[j CN 
relations by comparing the critical density of the gas tracer to the median 
density of the ISM. 



in dense star-forming environments. The circles and crosses 
correspond to the 7=1-0 and the 7=2-1 rotational transitions 
of HCN, respectively. 

The Lir-L' hcn relation (Figure O is tight and linear over 
several orders of magnitude, implying that stars form out of 
dense gas (Gao & Solomon 2004a,b; Wu et al. 2005). The 
dense gas fraction (L' HCN /L' C0 ) is nearly constant for galax- 
ies with LpR < 10 n L Q (Gao & Solomon 2004b), so Lqq can 
be used to indirectly trace the dense gas mass M* nse . How- 
ever, the dense gas fraction increases dramatically in LIRGs 
and ULIRGs (Lir > 10 u Lq), so CO does not trace dense gas 
mass in these galaxies (Gao & Solomon 2004b). HCN, on the 
other hand, has a critical density for excitation that is ^2 or- 
ders of magnitude larger than that of CO, so it traces dense, 
optically thick gas in star-forming GMC cores rather than dif- 
fuse GMC envelopes. The dynamical time for HCN-emitting 
gas is much less than the main-sequence lifetime of the most 
massive stars, so the intermittency factor for HCN-emitting 
gas will be approximately unity: 

<g™ ~ 2 x 10 5 yr ^ t /2 HCN < t MS -> £ « 1, (13) 

where p c rit,HCN ~ 10" 19 g cm" 3 . 

If the picture of radiation pressure feedback is correct, then 
it should determine the Lir-L hcn correlation directly. In fact, 
both the Eddington limit and the data show a linear relation 
between Lir and L HCN . The galaxies closely follow but do 
not exceed the Eddington limit for our preferred value of the 
Rosseland-mean opacity (kftr = 5-10 cm 2 g" 1 /dg.iso)- If the 
opacity is higher (kfir = 30 cm 2 g" 1 /dg.50), then many galax- 
ies are consistent with Eddington and a number are super- 
Eddington. For any of the values of the opacity that we as- 
sume, the general agreement between Lir and L' HCN suggests 
that radiation pressure may play an important role in regulat- 
ing star formation (Scoville et al. 2003). However, a number 
of important factors remain uncertain, which we discuss in $4] 

3.2. Molecular Schmidt Law and Radiation Pressure 
The Schmidt law is a tight power law relation between the 
surface density of star formation rate E* and the gas surface 
density (E* oc E^ 4 ; Kennicutt 1998). Furthermore, Bigiel et 
al. (2008) found that the Schmidt law for molecular gas is lin- 
ear within local star-forming galaxies (E* oc E H °). In the left 

panel of Figure[3] we plot E+ vs. Eh 2 for individual apertures 
of THINGS galaxies (small dots; Bigiel et al. 2008; Leroy 
et al. 2008), THINGS galaxies with H2 detections (open cir- 
cles), starburst galaxies 1 1 (solid circles), M82 super star clus- 
ters (stars; McCrady et al. 2003; McCrady & Graham 2007), 
and the Galactic Center star cluster (diamond; Paumard et 
al. 2006). We compare the data with the Eddington limit us- 
ing E* and Eh, as proxies for the radiation and gravitational 
forces, 



ecK 

11 Aravena et al. 2008; Becklin et al. 1980; Benford et al. 1999; Capak 
et al. 2008; Casoli et al. 1989; Chapman et al. 2005; Coppin et al. 2009; 
Coppin et al. 2010; Daddi et al. 2009; Downes & Eckart 2007; Downes & 
Solomon 1998; Greve et al. 2005; Knudsen et al. 2007; Mauersberger et 
al. 1996; Momjian et al. 2007; Neri et al. 2003; Paumard et al. 2006; Riechers 
et al. 2007; Riechers et al. 2008; Sajina et al. 2008; Sakamoto et al. 2008; 
Schinnerer et al. 2006; Schinnerer et al. 2007; Schinnerer et al. 2008; Smith 
& Harvey 1996; Walter et al. 2003; Walter et al. 2009; WeiB et al. 2001; Yan 
et al. 2010; Young & Scoville 1982; Yun et al. 2001 
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gas. The dot-dashed lines (both panels) are the intermittent Eddington limit (eq.[8j. The hatched regions (both panels) are the critical surface density or pressure 
for h = 30-100 pc (eq.QJj} where r M s ~ 2« dyn and rp IR ~ 1. S H , was calculated fromL^ Q using if L IR < 10 11 L andXj™ if Lir > 10 11 Lq. Overall, 
the Eddington limit suggests that radiation pressure sets an upper bound to the E» or the radiation pressure of a star-forming region or galaxy. Most star-forming 
regions or galaxies are sub-Eddington, but a few THINGS apertures and optically thick starbursts are super-Eddington (for k = 10 cm 2 g _1 /dgjso)- Several more 
optically thick starbursts will be consistent with or even exceed the Eddington limit for k = 30 cm 2 g _1 /dg.50- The rough agreement between starburst galaxies 
and the intermittent Eddington limit reinforces the likely importance of intermittency. However, the intermittent Eddington limit mildly under-predicts and 
P ra d at for Eh 2 < 10 Mqpc~ 2 and P m jd < 10 -12 ergs cm~ 3 , indicating that the effect of intermittency may be overestimated. 



where e«5x 10~ 4 is the efficiency of the conversion of mass 
into luminosity during the star formation process assuming 
a Kroupa (2001) broken power law IMF that extends up to 
120 M Q . For star clusters we assume a light-to-mass ratio 
appropriate for a zero age main sequence stellar population 
(vf = 3000 ergs s _1 g _1 ) and that the stellar mass is a lower limit 
on the gas mass of the parent GMC. We use the same Xqo 
values as in Figure [Hand again caution that Xco is uncertain 
to a factor of a few and may vary systematically from normal 
galaxies to starbursts (see 34.51 ). 

For the molecular Schmidt law, we find that the Edding- 
ton limit is an upper bound to E*. Most star-forming regions 
and galaxies follow the Eddington limit (solid line) and are 
within ~1.5 dex of it. The Eddington limit accounting for 
intermittency (dot-dashed line) appears to agree better with 
the data than the naive single-scattering Eddington limit, sug- 
gesting that intermittency may be an important effect. As 
the medium becomes optically thick near the critical surface 
density for intermittency (hatched region; see 32.21 ). the op- 
tically thick Eddington limit (E* dd oc E H2 /kfir oc E^°) pro- 
vides a firm upper bound to E* for our preferred value of the 
dust opacity (kfir = 5-10 cm 2 g" 1 /d g ,iso). If the dust opacity 
(^fir = 30 cm 2 g" 1 /d g ,5o) is enhanced due to a higher assumed 
dust-to-gas ratio, then some galaxies reach the optically thick 
Eddington limit and a few galaxies exceed it. 

In the right panel of Figure [3] we plot the radiation pres- 



sure from UV and FIR photons versus the midplane pressure. 
These pressures will balance each other at Eddington (solid 
line): 

F 7T 
Pmd ~ (1+TfTr)— ~ P m id = -GEgStot, (15) 

c 2 

where we take E tot = E g + 0.1E+ (Wong & Blitz 2002). The 
^rad-fmid plot shows that radiation pressure correlates strongly 
with midplane pressure over 10 orders of magnitude. The 
Eddington limit serves as a rough upper limit to f ra d, and 
most galaxies are within 2 dex of the Eddington limit. We 
note that some of the THINGS apertures and some dense star- 
bursts reach or exceed the Eddington limit. For galaxies with 
Pmid < P g-it, the critical midplane pressure (hatched region; 
see 32.21 & 14.21 ). we expect that the effects of intermittency 
are important; however, the intermittency adjusted Eddington 
limit (dot-dashed line) under-predicts P ra d for star-forming re- 
gions with P m id < 10~ 1L5 ergs cm" 3 . The intermittency factor 
may overestimate the importance of intermittency because of 
the si mpli fying assumption that subregions are "on" or "off" 
(see 132.21 ). We also see that galaxies and star-forming re- 
gions with 10 -11 ergs cm" 3 < P m jd < P C nt tend to fall signifi- 
cantly below the Eddington limit, possibly because our simple 
parametrization of Xco (see 34.51 ) is overestimating Mh 2 (and 
^mid) for these systems. Radiation pressure becomes increas- 
ingly more important in the optically thick limit (P ra d > P a it) 
as some galaxies and star-forming regions meet and exceed 
Eddington. As expected from Figures [T] & [2] if we assume 
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a larger dust-to-gas ratio and opacity (k = 30 cm 2 g" 1 /d g ,5o) 
potentially appropriate for dusty galaxies, then more of the 
optically thick starbursts would be super-Eddington. 

3.3. Radiation Pressure on Sub-galactic Scales 

So far we have evaluated radiation pressure on a galaxy- 
wide scale; however, the distribution of gas and star forma- 
tion in galaxies is inhomogeneous. Consequently, the Ed- 
dington ratio (r = frad/^mid) will likely vary on sub-galactic 
scales. We use observations from the THINGS survey (Wal- 
ter et al. 2008; Leroy et al. 2008; Bigiel et al. 2008) to calcu- 



late the Eddington ratio as a function of radius in azimuthally 
averaged radial bins and for semi-resolved (750 pc) aper- 
tures. Since the THINGS galaxies are generally in the single- 
scattering limit (see eq. |3), we conservatively adopt the radia- 
tion pressure to be P^ = Fir/c (see eq.fTBTl. For the midplane 
pressure given in eq. [O] the corresponding Eddington ratio 
is r tot = P^ d /(Q. 57rG£g£tot). Stars and atomic gas may not 
contribute significantly to the surface density in regions of ac- 
tive star formation, so we also calculate the Eddington ratio 
assuming that the midplane pressure depends only on the to- 
tal gas surface density T g = P^ d / {nGYfy or the molecular gas 
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FIG. 5. — The molecular gas Eddington ratio Th 2 = / 3 rad/(i"GSf Ii ) as a 
function of radius for NGC 6946. The line styles show the uncorrected Ed- 
dington ratio (triangles and thin solid line) and Th 2 corrected for an Xco 
gradient (dashed li ne), for a dust-to-gas ratio gradient plus a factor of <f> 2 
(<fi — h/RouC, see 32.21 dotted line), and for Xco and dust-to-gas ratio gradi- 
ents plus the intermittency factor and the cfr factor (thick solid line). After all 
of these factors are accounted for, the Eddington ratio is ~1 and nearly flat 
as a function of radius. We only show the effects of the Xco and dust-to-gas 
ratio gradients on the Eddington ratio in NGC 6946, but the profiles from the 
other THINGS galaxies shown in Figure|4]are qualitatively similar. 

surface density Tjj 2 = P^/ (wGT^ ). Intermittency may be 
important because the THINGS observations cannot resolve 
individual star-forming regions. We calculate the Eddington 
ratio corrected for intermittency (r{"{ = r tot /£, see eq. [HJ. In 
Figure|U we plot r tot (open circles), Th, (solid triangles), and 
rjJJJ (open squares) as a function of radius for azimuthally av- 
eraged rings. We find that T g is similar to r tot , so we omit T g 
for clarity. 

At intermediate radii (r ~ 1 — > several kpc), r tot and Th, 
generally increase from sub-Eddington (r ~ 0. 1 ) to approach- 
ing or exceeding the Eddington limit (r ~ 1). r tot reaches a 
maximum Eddington ratio at r ~ 5-10 kpc, where it falls off 
steeply. As the observations near the H2 detection threshold, 
Th, increases rapidly due to a small £h 2 with large error bars 
(see, e.g., Figure 40 from Leroy et al. 2008), and thus the large 
value of Th, at large r is consistent with Eddington to within 
the errors on £h 2 ■ For r > 1 kpc where r tot < 1 , the intermit- 
tency factor can boost rj£{ to the Eddington limit, suggesting 
that intermittency is important. 

The r < 1 regions in the inner parts of star-forming galax- 
ies present a challenge for radiation pressure regulated star 
formation. Intermittency cannot account for the low Edding- 
ton ratios of these regions. However, a metallicity gradient, as 
seen in observations of star-forming galaxies (Munoz-Mateos 
et al. 2009), increases the Eddington ratio at small radii (see 
jj4.41 i. A metallicity gradient that rises at smaller radii corre- 
lates with a decreasing gradient in Xco (Sodroski et al. 1995; 
Arimoto et al. 1996) and an increasing gradient in the dust-to- 
gas ratio (Munoz-Mateos et al. 2009). We adopt the Xco gra- 
dient given by eq. 10 of Arimoto et al. (1996) for data from the 



FIG. 6. — Histogram of the Eddington ratio (V = irad/flnid) f° r 750 pc 
apertures of THINGS galaxies. The midplane pressure was calculated us- 
ing two different methods: Th 2 = Prad/(7rG£g ) (dashed line) and r to t = 

P rad /(0.57rGSgE to t) (dotted line). The solid line (I**) shows the effect of 
adjusting the Eddington ratio by the intermittency factor (eq. [§J. r to i was 
calculated for all apertures with either an H2 or HI detection, but Th 2 could 
be calculated only for apertures with H2 detections. The r to t distribution is 
mostly sub-Eddington, but intermittency pushes the majority of the r]"j dis- 
tribution above the Eddington limit, implying that £ may overestimate the 
importance of intermittency. The Th, distribution is less peaked and shifted 
to higher Eddington ratios than the r to t distribution. Some apertures in the 
Th, distribution are at or above the Eddington limit in spite of the observa- 
tions not being able to resolve individual star-forming regions. 

Milky Way, M31, and M51 QogX /X e = 0.41[r/r e - 1], where 
r e is the effective radius, which we assume to be 7 kpc and 
X e is the value of Xco at the effective radius). To account for 
the dust-to-gas ratio gradient, we use a power law interpola- 
tion between /d g = 1/30 at 0. 1 kpc and /d g = 1/150 at 10 kpc, 
motivated by Figure 15 of Munoz-Mateos et al. (2009). In 
addition, collapsing GMCs e nhan ce the surface density by a 
factor of (f) 2 (4> = h /Rgmc\ see l2.2K making some regions opti- 
cally thick to FIR radiation. In Figure|5] we show the molec- 
ular Eddington ratio as a function of radius for NGC 6946 
since this galaxy is well below (~2 dex) the Eddington limit 
at small radii (see Figure [4j- After accounting for intermit- 
tency, an Xco gradient, a dust-to-gas ratio gradient, and a sur- 
face density enhancement in the GMCs, we find that Th, ~ 1 
for almost all radii in NGC 6946. We find qualitatively simi- 
lar results for all of the THINGS galaxies shown in Figure [4] 
assuming that the metallicity, Xco, an d dust-to-gas ratio gra- 
dients are similar to those adopted for NGC 6946. Thus, it 
is at least in principle possible to explain the nominally sub- 
Eddington inner regions of local star-forming galaxies using 
a combination of these effects. 

In Figure |6l we plot the distributions of the individual Ed- 
dington ratios for the THINGS apertures (750 pc resolution): 
r tot (dotted line), r H , (dashed line), and rj"| (solid line). The 
distribution of r tot is peaked around r tot ~ 0. 1 , and the major- 
ity of apertures are sub-Eddington for r tot . The high r tot tail 
of the distribution extends above the Eddington limit, with 
super- Eddington apertures comprising 5% of the total aper- 
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tures and containing 5% of the total flux. Star-forming regions 
are unresolved on 750 pc scales, so r tot should be adjusted 
to account for intermittency (T{"{). However, most apertures 
lie above the intermittency adjusted Eddington limit. As in 
Figures Q~|-|5] this shift suggests that intermittency is impor- 
tant for radiation pressure supported star formation in normal 
spirals; however, £ appears to overestimate the importance of 
intermittency, possibly due to the simplifying assumption that 
subregions are either "on" or "off" (see 32.2l & l4~2l i. The dis- 
tributions of r tot and r g are similar, so we do not plot T g for 
clarity. 

The distribution of Th, is less peaked and shifted to sys- 
tematically higher values than the distribution of r tot with 
10% of these apertures radiating at or above the Eddington 
limit. This is not surprising (given Figure because radia- 
tion pressure will likely be more important in H2-dominated 
star-forming regions. The detection limit for H2 is higher than 
that for HI, so the distribution of Th, contains fewer aper- 
tures than r tot . In addition, the apertures with H2 detections 
tend to be within 0AR 2 s (Bigiel et al. 2008), so they might 
have increased metallicities and dust-to-gas ratios with de- 
pressed Xco values, which would increase the Eddington ratio 
(see Figure|5] jj4.4l and l4.51 l. Super-Eddington apertures con- 
tain 6% of the total flux in H2-detected apertures across the 
whole sample; in NGC 6946, for example, super-Eddington 
apertures contain 10% of the total flux. The super-Eddington 
apertures indicate that radiation pressure can be dynamically 
dominant even when individual star-forming regions remain 
unresolved, suggesting that radiation pressure may be more 
important on the scale of GMCs and massive star clusters. Fi- 
nally, we calculated Th, assuming that UV photons contribute 
to the radiation pressure f rac j = (Fuv+Fir)/c. The distribution 
of Fh 2 remains nearly the same because star-forming regions 
have Fir/Fuv >1,so we refrain from plotting it in Figure [6] 

4. DISCUSSION 

We have compared globally-averaged and resolved obser- 
vations of star-forming galaxies with theoretical expectations 
based on the theory of radiation pressure supported star 
formation (see §2). Although the uncertainties are large (see 
below), our primary findings are as follows. 

1. Figures QJO show that star-forming galaxies meet, but 
do not dramatically exceed, nominal expectations for the 
dust Eddington limit. When some subregions do seem to 
exceed the Eddington limit (as in the outer regions of galaxies 
in Figure [4] & [5), we consider this to be consistent with 
Eddington since trends in the dust-to-gas ratio and CO-to-H2 
conversion factor, as well as the large-scale stellar potential 
and intermittency of the star-formation process (£; eq. [8} 
affect the Eddington ratio at order unity. 

2. The Lir-L' hcn plot (Figure I3 provides the strongest 
evidence for the importance of radiation pressure feedback 
since L' HCN is expected to directly trace the dense actively 
star-forming gas and Lir traces the total star formation rate. 
If radiation pressure in fact dominates feedback, we would 
expect a one-to-one correspondence between these two 
quantities, and such a relation is observed (see also Scoville 
et al. 2001; Scoville 2003). Nevertheless, for typical values 
of both k in the optically thick limit and the HCN-to-H2 
conversion factor, the Eddington limit overpredicts Lir by 
a factor of ~3-6. This discrepancy may indicate that the 
dust-to-gas ratio is larger in the dense HCN-emitting regions, 



or that the HCN-to-H2 conversion factor is smaller (see 
34.51 below). If radiation pressure feedback regulates star 
formation, then this relation is in a sense more fundamental 
than the Schmidt Law because HCN-emitting gas is more 
closely connected with star formation than CO-emitting gas, 
in which case we would expect S^ dd = 47rG£ d i ense /(ecKFiR). 

3. The central regions of all galaxies in Figure |4] are prima 
facie substantially sub-Eddington when a constant dust-to-gas 
ratio and CO-to-H2 conversion factor are applied to all sub- 
regions without regard to their radial location. If radiation 
pressure is in fact the dominant feedback mechanism in these 
regions, a much higher central dust-to-gas ratio and a lower 
CO-to-H2 conversion factor are required (see Figures [4]&[5]l. 
It would be particularly useful for testing radiation pressure 
feedback to produce the same profiles in HCN. 

4. The "break" in the observed Schmidt Law at £ g ~ 100- 
1000 Mq pc -2 (see Figure [3j Daddi et al. 2010b) may be 
due to the transition from the single-scattering limit to the 
optically thick limit in the GMCs that collapse to form stars, 
as in MOT. 

5. If radiation pressure is the primary feedback mechanism 
for regulating star formation, then we predict that the Schmidt 
Law will follow the form of eq. fT4l(for d iscussion of k and £ 
as well as uncertainties see ^2l& l4.2144.5b . 

6. A testable prediction of radiation pressure feedback is that 
all else being equal the star formation rate should depend 
linearly on the dust-to-gas ratio in the optically thick limit. 

In addition to these points, below we note an implication of 
radiation pressure feedback that has so far not been stated in 
the literature ( 134. 11 1. Finally, in the remaining subsections we 
highlight the dominant uncertainties in our work as a guide 
for future research on the importance of radiation pressure 
feedback in star-forming galaxies. 

4.1. Gas Depletion Timescale 

The gas depletion timescale, the time required to consume 
a galaxy's gas reservoir at the current SFR, is observed to be 
~2 Gyr in normal spirals (Kennicutt 1998; Leroy et al. 2008). 
Radiation pressure sets the gas depletion timescale to be 

_ M g _ M g ec 2 _ ecu 

'gas = T^ - — "77 — ~j 7-T7- (10) 

M+ £L E dd 4ttG£ 

Using typical numbers for a spiral galaxy in the single- 
scattering limit, the gas depletion timescale is 

f gas ~ 2.1 Gyr 5^ /2 /^i;, (17) 

where E30 = S/30 M Q pc" 2 and h\oo - h/100 pc. This nor- 
malization of ? gas is in good agreement with the observed gas 
depletion timescale, but eq. [TTJpredicts that the gas depletion 
timescale should have a strong dependence on S g , in contrast 
to the observations of Leroy et al. (2008) (see their Figure 
15). For completeness, we note that variations in the dust-to- 
gas ratio will not affect the dependence of f gas on £ g in the 
single-scattering limit, but uncertainties in Xco, £, and <\> (see 
34.4144. 5b might impact the gas depletion timescale. 

Hot starbursts and optically thick subregions (see 32.1.2b 
have intermittency factors that approach unity and nearly con- 
stant opacities, so the gas depletion time is approximately 



10 



constant, 

f gas «5.7MyrK 10 , (18) 

where kio = kfir/10 cm 2 g" 1 . For comparison, Sakamoto et 
al. (2008) find that the optically thick western nucleus of Arp 
220 has a gas depletion time of —6 Myr. The fact that the gas 
depletion timescale set by radiation pressure is consistent with 
the observed gas depletion timescale in spirals and ULIRGs 
is equivalent to the statement of Figures QJ[5] that starbursts 
approach the dust Eddington limit. In addition, we point out 
that radiation pressure feedback predicts that the specific SFR 
(SSFR) will be SSFR - 47rG£/ gas /(ec>t) for small / gas . 

4.2. Intermittency 

The intermittency factor £ (see 32.21 MQT) relates the 
properties of radiation pressure dominated star-forming sub- 
regions to the global properties of a galaxy. However, £ may 
overestimate the effect of intermittency in some galaxies (see 
Figures [3] [4] & [3}. We expect the determination of £ to be 
complicated by uncertainty in the timescale for the central 
cluster of a sub-unit to be bright (/ms ~ 4 Myr). We adopt 
?ms as the time a cluster will be bright, since the cluster lumi- 
nosity drops rapidly after the most massive stars in the cluster 
explode as supernovae. However, a cluster will continue to 
emit after ?ms- Indeed, models of cluster luminosity indicate 
that a similar amount of momentum will transferred to the 
gas during the time to — > f M s and during the time f M s — > 4f M s, 
where the cluster luminosity has dropped by 1 dex after ~4?ms 
(Leithereret al. 1999). Further uncertainty in £ is due to ambi- 
guities in calculating the lifetime of a sub-unit (— 2fd yn + fMs), 
especially the disk dynamical time (see eq. |9). The dynami- 
cal time likely varies with galactocentric radius because £ is 
a strong function of radius (Leroy et al. 2008). Thus, trying 
to determine an effective £ applicable to a galaxy as a whole 
may be difficult if £ changes locally. Overall, we expect the 
uncertainty in £ to be a factor of a few to several. 

4.3. The FIR Optical Depth 

A key theoretical uncertainty in calculating the Eddington 
limit for dense starbursts is the effective optical depth (r e ff) for 
surface densities where tfir > 1 . In order for radiation pres- 
sure to be dynamically important in optically thick GMCs, 
r e ff must exceed unity. Based on the high Mach number tur- 
bulence simulations of Ostriker et al. (2001), MQT conclude 
that if the ISM is optically thick on average, then the vast 
majority of sight lines will be optically thick. For compari- 
son, KM09 argue that instabilities, such as Rayleigh-Taylor 
and photon-bubble instabilities, will reduce the effective opti- 
cal depth of the dense ISM to ~ 1 . However, MQT note that 
both the midplane pressure from gravity P m jd ~ ttGY, 2 and op- 
tically thick radiation pressure P lat j ~ tF jc oc Y? scale as E 2 , 
a feature unique to radiation pressure among stellar feedback 
processes. Thus, if radiation pressure cannot regulate star for- 
mation in dense, optically thick gas, then no known stellar 
feedback process can. 

4.4. Dust-to-Gas Ratio and Metallicity 

The coupling between radiation and gas directly depends 
on the dust-to-gas ratio (k oc /d g ). In this analysis, we as- 
sume the Galactic value for the dust-to-gas ratio (/d g = 1 /150) 
and solar metallicity; however, there is strong evidence that 
/d e and metallicity change with environment. The dust-to- 
gas ratio has been shown to correlate with metallicity and 
radius (Issa et al. 1990; Lisenfeld & Ferrara 1998; Draine 



et al. 2007; Munoz-Mateos et al. 2009). Munoz-Mateos et 
al. (2009) find that the dust-to-gas ratio can climb to values as 
high as /d g ~ 1/10 in the centers of spiral galaxies. This in- 
crease in metallicity and dust-to-gas ratio is necessary for the 
centers of star-forming spirals to be at Eddington (see Figures 
H&Qand 33.31 l. The average dust-to-gas ratio of local spirals 
also varies by a factor of a few (e.g., M51 has a /dg ~ 1 /75; 
Draine et al. 2007). Furthermore, the dust-to-gas ratio is ob- 
served to be higher in some dense starbursts, such as SMGs 
(/d g ~ 1/50; Kovacs et al. 2006; Michalowski et al. 2010) 
and sub-mm faint ULIRGs (/ dg ~ 1 /20; Casey et al. 2009). 
Importantly, if we adopt a dust-to-gas ratio potentially appro- 
priate for dusty starbursts (short dashed line in Figure Q] [2] 
and the left panel of Figure |3j, then a substantial fraction of 
optically thick galaxies would be at or above the Eddington 
limit (Figure|5]illustrates this for NGC 6946). 

4.5. Molecular Gas Tracers 

The Eddington limit depends strongly on the CO-to-H2 
(Xco) and the HCN-to-H2 (^hcn) conversion factors (see eqs. 
fm & H2l . These conversion factors are two of the largest 
sources of observational uncertainty in our calculations be- 
cause they vary with excitation conditions (X oc ^np^/Tb, 
where Tf, is the brightness temperature) and metallicity. Sev- 
eral lines of evidence suggest that X^q overestimates Mh, in 
starburst galaxies. For example, Xco is a factor of —3 lower 
in M82 (Weiss et al. 2001) and a factor of —5 lower in a 
sample of local ULIRGs (Downes & Solomon 1998). To ac- 
count for the different values of Xco in normal spirals and ex- 
treme starbursts, we apply the Milky Way Xco value to galax- 
ies with Lir < 10 n L© and the ULIRG Xco value to galaxies 
with Lir > 1O H L0. Because this prescription is somewhat 
simplistic, it probably overestimates Mh, in moderate lumi- 
nosity starbursts (L IR < 10 u Lq), such as M82, and in the 
centers of star-forming spirals (see Figures [4] & [5] and 33.31 ). 
Additionally, it likely underestimates Mh 2 in ultra-luminous 
(Lir ~ 10 12 L0) high redshift disk (BzK) galaxies, for which 
Daddi et al. (2010a) find a value of Xco that is consistent with 
the Galactic value. As a result, moderate luminosity starbursts 
and the centers of star-forming spirals may be closer to Ed- 
dington and BzK galaxies might be further below the optically 
thick Eddington limit than Figure [3] would suggest. 

Unfortunately, Xhcn is more uncertain than Xco because 
there is no direct calibration of Xhcn from Milky Way GMCs. 
For normal spirals, Gao & Solomon (2004a,b) find Xhcn ~ 
10 M Q (K km s" 1 pc 2 )" 1 for virialized cloud cores with (n) = 
3 x 10 4 cm" 3 and 7J, = 35 K. They caution that X H cn could be 
lower in regions of massive star formation due to significantly 
higher brightness temperatures 7b ~ few x 10 2 K (Boonman 
et al. 2001). Ultra-luminous starbursts exhibit widespread in- 
tense massive star formation, so one might expect that Xhcn 
is lower in more luminous galaxies. For example, Gracia- 
Carpio et al. (2008) estimate that Xhcn should be —4.5 times 
lower for galaxies at Lfir ~ 10 12 L© than at Lfir ~ 10 11 L Q . 
We note that if Xhcn is smaller than the assumed value of 
3 M Q (K km s" 1 pc 2 )" 1 , then more galaxies will approach or 
exceed the Eddington limit (see Figure 13. For example, de- 
creasing Xhcn by a factor of — 2 brings essentially all galaxies 
in line with the Eddington limit for the nominal value of kfir 
(-5-10). 
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